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Abstract 

The two-dimensional scaling Ising model in a magnetic field at critical temperature is in- 
tegrable and possesses eight stable particles Ai [i = 1,...,8) with different masses. The 
heaviest five lie above threshold and owe their stability to integrability. We use form factor 
perturbation theory to compute the decay widths of the first two particles above threshold 
when integrability is broken by a small deviation from the critical temperature. The lifetime 
ratio ^4/^5 is found to be 0.233; the particle A5 decays at 47% in the channel AiAi and for 
the remaining fraction in the channel A1A2. The increase of the lifetime with the mass, a 
feature which can be expected in two dimensions from phase space considerations, is in this 
model further enhanced by the dynamics. 



1 Spectrum and analyticity in Ising field theory 



Quantum field theory in 1 + 1 dimensions allows for a subclass of integrable cases. While inte- 
grability is hardly detectable at the level of correlation functions, it has very visible implications 
at the scattering level. Indeed, the presence of infinitely many conserved quantities forces the S- 
matrix to be completely elastic and factorized pp. These properties, in turn, allow for a feature 
of the particle spectrum which is peculiar to integrable quantum field theories: the existence of 
particles with mass higher than the lowest threshold which are stable even in the absence of an 
internal symmetry which would prevent them from decaying. Once integrability is broken by 
an arbitrarily small perturbation, the stability of these particles would contradict though the 
analyticity requirements on the ^-matrix, so that they are forced to develop a finite lifetime. 
This phenomenon already appears in the field theory describing the simplest universality class 
of critical behavior, namely that of the Ising model in a magnetic field [Ij. The quantitative 
study of the decay processes following the breaking of integrability in Ising field theory is the 
subject of this paper. 

The scaling Ising model in two dimensions is described by the action 

A = Ao-T j (fx£{x)-h j (fxa{x), (1) 

in which the critical point action is perturbed by the energy operator e{x), with scaling 
dimension = 1, and by the spin operator a{x), with scaling dimension X^- = 1/8. The 
coupling constants^ 

r ~ M^-^- = M , 
h ~ M^-^- = M^^/^ 

measure the deviation from critical temperature and the magnetic field, respectively, and enter 
the dimensionless parameter 

^ = ^ (2) 

which can be used to label the renormalization group trajectories flowing out of the critical 
point. 

It is well known that the action with h = describes free excitations with fermionic 
statistics. These are ordinary particles in the high-temperature phase r > (ry = +00), and 
topological excitations (kinks) interpolating between the two degenerate vacua of the r < 
phase {rj = —00) in which the spin reversal simmetry is spontaneously broken. The evolution 
of the particle spectrum at /i 7^ was first discussed in [21. In the limit rj —00 the particle 
spectrum exhibits infinitely many discrete levels generated by the confinement of the kink- 
antikink pairs induced by an arbitrarily small magnetic field removing the degeneracy among 
the two vacua of the low-temperature phase. Hence, the spectrum in this limit is reproduced by 



^The mass scale M is the inverse of the correlation length. 
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the non-relativistic result for a linear confining potential^ 

mn ~ 2m H -^^ — , (3) 

where m is the mass of the kinks, v the spontaneous magnetization and Zn, n = 1,2,..., are 
positive numbers determined by the zeroes of the Airy function, Ai(— z„) = 0. Of course, the 
particles with masses m„ larger than twice the lightest mass mi are unstable. It was argued 
in |5] that the number of stable particles descreases as rj is increased, until a single asymptotic 
particle is left in the spectrum at r; = +oo. This general pattern has been confirmed by numerical 
investigation of the spectrum of the field theory [1JI5]. In particular, it has been found that 
the second and third lightest particles become unstable at 

r?2 = 0.333(7) , % ^ 0.022 , (4) 

respectively These values of r/, as well as the other normalization-dependent numbers we 
quote in the following, refer to the normalization of the operators in which 

(a(x)(j(0)) ^ 

(e(x)e(O)) - (5) 

as \x\ 0. Corrections to the mass mi for t] +00 as well as to the spectrum (jSJ have also 
been computed IH] • 

A. Zamolodchikov discovered at the end of the eighties that the Ising field theory is 
integrable for r = and computed the exact S'-matrix through the bootstrap method j2- He 
found, in particular, that the spectrum of the theory on the magnetic axis contains eight stable 
particles Aa (a = 1, . . . , 8) with masses^ 

mi = (4.40490857..) \hf^^^ 

m2 = 2mi cos - = (1.6180339887..) mi 
5 

ms = 2mi cos— = (1.9890437907..) mi 
30 

m4 = 2m2 cos— = (2.4048671724..) mi 
2vr 

m5 = 2m2Cos — = (2.9562952015..) mi 
15 

me = 2m2 cos = (3.2183404585..) mi 
ou 

vr 77r 

m^ = 4m2 cos — cos — = (3.8911568233..) mi 
5 30 

mg = 4m2 cos — cos — = (4.7833861168..) mi . 
5 15 



^The spectrum @ was originally obtained in through the study of the analytic structure in momentum 
space of the spin-spin correlation function for small magnetic field. 

^The constant relating mi to the magnetic field was determined in |7|. 
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Notice that the last five masses all lie above the lowest two-particle threshold 2mi . Since the Ising 
model in a magnetic field possesses no internal symmetries (so that all particles are completely 
neutral), the decay processes Aa — > AiAi, a = 4,... ,8, are allowed both by kinematics and 
by symmetry. The stability of the particles above threshold, however, does not violate any 
requirement of the S'-matrix theory as long as integrability is present. 

To see this, consider for example the scattering amplitudes of Ai with itself and with A2. 
They read"^ |2] 

Su{e) = h/3{e)t2/,ie)h/i5io) (6) 

Si2{0) = i4/5(^)%5(^)*7/15(^)Vl5 (0) , (7) 

with 

_ tanh^{e + ina) 

to^y^l - - — T-TTa '■ — ^ ■ 

tanh 2\" ~ t-T^oi) 

The simple poles at = iu^^ in the amplitude Sab{0) indicate that the particle A^. with mass 
square 

ml = + ml + 2mami, cos (9) 

appears as a bound state in the scattering channel AaAi,. Hence we see that the scattering chan- 
nel AiAi produces the first three particles as bound states, while the channel A1A2 produces the 
first four. Figures l.a and l.b show the corresponding analytic structure for the two amplitudes 
in the complex plane of the Mandelstam variable s (square of the center of mass energy). As 
required by the S'-matrix theory, all the poles corresponding to stable bound states lie on the 
real axis below the lowest unitarity threshold in the given scattering channel. This is true for 
all the amplitudes Sab, a, 6 = 1, . . . , 8, of the Ising field theory at r = (see e.g. |H] for the full 
list of amplitudes). When integrability is broken (i.e. as soon as we move away from r = 0), 
however, the inelastic channels and the associated unitarity cuts open up. In particular, the 
process A1A2 — > ^lAi acquires a non-zero amplitude, so that the threshold located at s = 4m^ 
becomes the lowest one also in the A1A2 scattering channel. Since the pole associated to A4 is 
located above this threshold, it can no longer remain on the real axis, which in that region is 
now occupied by the new cut. The position of the pole then develops an imaginary part which, 
according to the general requirements for unstable particles [2], is negative and brings the pole 
through the cut onto the unphysical region of the Riemann surface (Fig. l.c). 



2 Mass corrections and decay widths 

The change in the position of the poles when the temperature is moved away from its critical 
value can be computed perturbatively in r within the framework of Form Factor Perturbation 

"^The energy and momentum of the particles are parameterized in terms of rapidities as {p",p^) = 
(rria cosh 6, TTia sinh 0). Relativistic invariant quantities as the scattering amphtudes depend on rapidity differ- 
ences only. 
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Figure 1: Poles and unitarity cuts for the scattering amplitudes and in the integrable 
case r = 0, (a) and (b), respectively, and for r slightly different from zero (c). In (c) the particle 
^4 became unstable and the associated pole moved through the cut into the unphysical region. 
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(b) 

Figure 2: Diagrams determining the leading corrections to the real (a) and imaginary (b) parts of 
the masses at small r. In (b) also the intermediate particles are on shell and energy-momentum 
is conserved at each vertex. For c > 5 also diagrams with more than two particles in the 
intermediate state contribute to the imaginary part. 



Theory around integrable models 1^. The form factors of an operator $(a;) are its matrix 
elements between the vacuum and the multiparticle asymptotic states. We denote them as 



FL.a„ (01, . . . , ^n) = {0\m\Aa, (^l) • • • (^n)) • 



(10) 



Then the leading corrections to the real and imaginary parts of the particle masses can be 
written as (see Fig. 2) 



5Rem2 ~2t/c, 



(11) 



Imm^ = — rricTc 

a<b , ma+mi,<mc 



E 



I fcab I 



a<b , ma+mi,<m,c ^c^a 

where Tc^ab is the decay width of the particle Ac into the channel Aa A^, and 



sinh^i^'^') 



(12) 



/,= F,^,(m,0)| 



r=0 



(13) 



/ cab 



cab 



^,^^e^cab)^Qicaby 



(14) 



The rapidities 9^'^^'^ and Q'^"' are determined by energy-momentum conservation for the decay 
Ac AaA{, in the rest frame of Ac- All the masses appearing in the r.h.s. of (|12|) are taken 
at r = 0. For c > 5 the sum in ()12() must be completed including the contributions of the 
decay channels with inorG tlictn two psirticlGs in the final state. Once the decay widths ^.q^^ 



^{cab) 
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are known, one can determine the lifetime tc of the unstable particle Ac, given by 

tc = ^ , = ^ Tc^ab • (15) 

a<b 

The form factor problem for the Ising field theory in a magnetic field at r = has been studied 
in jlUl lllj . where one and two-particle form factors have been computed, and these results are 
reviewed in [S]. We have now extended the form factor bootstrap to the three-particle case in 
order to compute the decay widths (|12jl . In the next section we recall the main steps of this 
approach and refer the reader to |HI for details. 



3 Form Factors 

Form factors in an integrable quantum field theory satisfy a set of equations |121 113j which in 
the case of neutral particles read 

^ai...a,a,+ i...a„ (6*1 , • • • , 6*1 , ... ,9n) = 

Sa,ai+iiOi - di+l)Fa^...ai+ia,...a„i^l^ ■ ■ ■ ■ ■ ■ , On) (16) 

F^,...a„iGl + 2^^, ^2, . . . , 0n) = i^*...a„ai (^2, • • • , ^n, ^l) (17) 

R'eS6i^-6»6=m=j F^bai...a„iOa,Ob,Oi, ... ,9n) = iglb Fcai...a„iOc, ...,0n) (18) 

Rese'=e+iTT Fabai...ani^' ^^^^^^ ■ • ■ ,6*71) = 

iSab (^l-U.Sa^-^^j-e)j F*...„J^i,...,^„) , (19) 

with the three-particle couplings g'^^^ entering ()18() determined from the scattering amplitudes 
through the relations 

Sab{9 Kb) ^ ■ (20) 

^ - Kb 

A parameterization of the form factors which solves (|16|) and I17|) and contains the poles pre- 
scribed by ((TH|l and (fT^ is given by 

pmin/Q, Q ,\ 

Here ^ 

Frr{e)=(-ismh^-^ " n (^7/3oW)'^ , (22) 

rr. f /■°°dtcosh(a-i)t ^{iTT-e)t\ 

satisfying 

To^io) = -ta{e)Tc.i-e) (24) 



with 
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T^{e + 2iTT) = T^{-e) 



(25) 



r„(0) ~exp 1^1/2 , 1^1 ^oo. (26) 

The quantities 7 and p-^ in 1)22(1 coincide with those entering the expression 

Sam= n (V30(^))"^ (27) 

for the amphtudes {Sab{^) = (— 1)*^"*). The dynamical poles are inserted in ()21() through the 
factors 

Dam = n (^7/30(e))^^ (^l-7/30(^))'^ , (28) 

l&Gab 

where 



cos vra — cosh 6 

and 



Te.{e) = (29) 



2 cos^ ^ 



= 71 + 1 , j^ = n, if = 2n + 1 
= n , = n , it p.y = Zn . 

Multiple poles corresponding to > 1 are related to multiscattering processes jTl] and give rise 
to residue equations which generalize ()18j) |lflj . 

The operator dependence in (PT|) is contained in the entire functions ...0,^(6*1, • • • , dn)- For 
scalar operators <^ with scaling dimension the asymptotic bound JU] 

hm F*...,„(^i, ...,9n)^ exp(y$|0,|) (31) 

< ^ , (32) 

together with the invariance under 9j — s- 6j + 2i7r, implies that the ...^^ are rational functions 
of the variables Xj = e^^ . In the two-particle case, taking into account Lorentz invariance, the 
vanishing of the residue (fTH)) and the identity -Faia2 (^i) ^2) = -^0201 (^i, ^2), one can write 

with 

pLMM = E 4'a2 cosh'^(^i - O2) . (34) 

fe=0 

In the three-particle case we use the notation 



ai,a2,a3=a'^,a'2,a'^ 
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The coefficients in (jSU and H35() are real and restricted in number by the bound (|32() . In 
particular, defining Deg„.„^, through 

and 

Degi_+ = -Tifigaiaj - Dega^afe " ^ataj - ^a.a^ (37) 

Degi _ = Deg„^„^. + Deg„^„^ (38) 

{i,j,k a permutation of 1,2,3), we conclude 

a[ = -Deg,,_ (39) 
a': = -Deg,,+ (40) 

for both $ = 0", e in 1)35^ . Further constraints on ()35|) come from symmetry under permutations 
of identical particles 

{01,02,03} _ (^*,{a2,oi,03} /^^N 

and from Lorentz invariance, which for scalar operators implies 

3 

k=l l<i<j<3 

The form factor equations together with the asymptotic bound allow to show that the space 
of non-trivial relevant scalar operators is two-dimensional [S], as expected for the Ising model, 
but cannot resolve the linear combination between a and e. This is achieved imposing the 
asymptotic factorization condition 



«5 



which holds, in particular, for relevant scalar scaling operators in theories without internal 
symmetries |15j . Of the two form factor solutions obtained in this way, that corresponding to 
a is picked up exploiting the proportionality between this operator and the trace of the energy- 
momentum tensor at r = |1U1 lllj : the remaining solution must then correspond to e. Once 
the initial conditions of the form factor bootstrap for <I> = cr, e have been fixed in this way, 
the functions Qq^ ..^^ in (|7T|) for these operators can be uniquely determined using the residue 
equations. We give in a list of results for n = 2 (which extends that given in lllj) and 
for n = 3, 4. 
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4 Results 



From the above results we extract the following values to be used in and ()12|) 

/i = (-17.8933..)(e).=o (44) 

/2 = (-24.9467..)(e)r=o (45) 

/3 = (-53.6799..)(e)r=o (46) 

U = (-49.3206..)(e)r=o (47) 

|/4ii| = (36.73044..) I (e)|,=o (48) 

1/511 1 = (19.16275..) 1(e) 1,^0 (49) 

|/5i2| = (11.2183..) . (50) 

Within the normalization @ the vacuum expectation value appearing in these results takes the 
value [Ej 

{e)r=o = (2.00314..) \hf/^' . (51) 

The variation of the mass ratios 

Re ma 

Ta = , (52) 

mi 

close to the magnetic axis, is given by 

Sra = -^(rl-^)+0{r'). (53) 



mima V " fi 

It can be immediately checked that the values H44() - H47|) imply positive coefficients for r in this 
equation when a = 1, ... ,4, in agreement with the expectation that the distance of the stable 
(unstable) particles from the threshold 2mi decreases (increases) as rj increases. The mass 
variations of the first three particles close to the magnetic axis were measured in |3 El and 
fully agree with the predictions of the Form Factor Perturbation Theory. 

It is tempting to ignore the non-linear corrections to H53|) and obtain a first-order estimate 
r]n^ of the values r/,„ at which the n-th particle decays. These are determined by the condition 

rn{Vn) = 2, n>l, (54) 

and, in this way, one obtains 

(55) 
(56) 
(57) 

Comparison with (jlj shows that the first-order estimate captures the correct order of magnitude 
of the decay point for n = 2,3. Notice that, if for n = 3 the mass ratio changes of only 0.5% 
from r/ = to 773, the variation of m-s and mi with respect to mi|^=o in the linear approximation 



= 2, 


n > 1 




~ 0.2 




~ 0.01 




~ -0.2 . 
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^ 

Figure 3: The first few trajectories wliich divide the h-r plane into regions with a different 
number of stable particles. There are n stable particles in between the trajectories labeled by rjn 
and rjn+i- These trajectories densely fill the plane when the negative orizontal axis is approached 
(?7 -oo). 



is considerably larger. The trajectories corresponding to r/2, % and to the first-order estimate 
for ri4 are shown in Fig. 3. 

For the imaginary parts of the first two particles above threshold we obtain 



Imml ~ (-840.172..) ( iMl^^ = (-173.747..) (58) 

nil J 

Imm'i ~ (-240.918..) ( = (-49.8217..) . (59) 



The ratio of lifetimes 

lim — = lim '""^ ^"^""^ = 0.23326.. (60) 



td rriA Im m? 

.. — = lim 1 

+0 r-»o ms Im m| 



is universal. While can only decay into AiAi, has also the A1A2 channel available. The 
relevant branching ratios 

"c^ab = Tf 2I (61) 

|im m^l 

are 

lim 65^11 = 0.47364.., lim 65_i2 = 0.52635.. . (62) 

It appears from H60|) that the lifetime of is more than four times longer than that of A4, 
and this seems to contradict, somehow, the expectation inherited from accelerator physics that 
the lifetime of an unstable particle decreases as its mass increases. Notice, however, that the 
d-dimensional phase space for the decay Ac — > AaAj, is 



/ 



o~ ^ yP<^ +Pb- Pc) ~ , (63 

PI Pb 



10 



where p = \pa\ = \pb\ is taken in the center of mass frame. For fixed decay products, p increases 
with rric and, in d = 2, it joins the factor nic in the denominator, leading therefore to a sup- 
pression of the phase space. The results (jl5|) - (|50|) show that such a suppression of the decay 
width is further enhanced by the dynamics in a way that is not compensated by the opening of 
additional decay channels. 

In we also give the four-particle form factor F^^]^^|t-=o which determines the first-order 
variation to the elastic scattering amplitude 5ii through the formula 

SSuiei - 62) ~ ^-—^ lim FfuiiOl + ivr + ,5, 02 + ivr + ,5, 61,62) . (64) 

smn(f7i — U2) o— »o 

Here the limit is needed to deal with the singularities prescribed by H19() . 



5 Numerical studies 

Checks of the analytic results for the particle spectrum may come, in particular, from the 
numerical diagonalization of the Hamiltonian on a cylinder geometry, the results for the plane 
being recovered as the circumference R of the cylinder goes to infinity. For the stable particles 
of the Ising model in a magnetic field this kind of analysis has been performed both in the 
continuum [^ElE] and on the lattice jl8 | I2() | 1^. always confirming the analytic predictions. 
The numerical tests are more complicated for unstable particles. 

Qualitatively, the signature of unstable particles on the cylinder is very clear [^. At the 
integrable point, when the energy levels are plotted as a function of R, the line corresponding 
to a particle above threshold crosses infinitely many levels which belong to the continuum when 
R = 00 (Fig. 4). Once integrability is broken, this line "disappears" through a removal of level 
crossings and a reshaping of the lines associated to stable excitations (Fig. 5). Quantitatively, 
though, it is not so obvious how the decay width can be measured from such a spectrum. In 
absence of a direct measurement of this quantity, one could try however to estimate, using the 
exact matrix elements computed above, the energy splitting resulting from the removal of a level 
crossing taking place at a sufficiently large value R* of the cylinder circumference. 

The usual formula for the first order energy splitting of two states |1) and |2), which are 
degenerate in the unperturbed system, is given by 

= V(lii - V22)2 + 4|Fl2|2 , (65) 

where Vij = {i\V\j) are the matrix elements of the perturbing Hamiltonian V. When we put 
the theory on the cylinder and perturb around r = 0, we have 

V = -T dx e{x) . (66) 
Jo 

For R* large enough, the corrections coming from the finite volume dynamics can be neglected in 
the first approximation. Hence, the only R dependence comes from the change of normalization 
of the particle states when passing from the plane to the cylinder: 

{Ab{pt)\AaiPa)) = 27r6abplSipl - pI) SabSpipiplR. (67) 



11 



By taking into account the normalization {i\j) = 6ij of the states entering ()65() . the matrix 
elements to be used for the removal of the crossing between the levels corresponding to ^4 and 
Ai^i are then given by 

where 

/nil = lim FUier^ +^. + 5, -6^^ + ^. + 5, -e^^''\ef''^) 

O— >0 T = 

= (408.78..) {e)r=Q ■ (71) 

For the purpose of comparison with numerical data it can be useful to consider the universal 
ratio AE/5Egs{R*), where 

6EgsiR) = -Rt {e)r=o (72) 

is the first order variation of the ground state energy. 

Energy spectra on the cylinder can be obtained directly in the continuum through the trun- 
cated conformal space approach [231, ™ which the Hamiltonian is diagonalized numerically on 
a finite dimensional subspace of the conformal states relative to the critical point. This is the 
method that we have used to obtain the numerical spectrum. By including states up to the 5th 
level of the Verma module of the conformal families of the Ising model, our final Hamiltonian 
has been truncated up to 43 states. Even though this truncation of the Hilbert space has proved 
to be successful to measure mass ratio corrections and the change of the ground state energy, it 
fails however to measure with sufficient precision the energy splittings which take place above 
threshold at large values of R. This was not totally unexpected since it is well known that 
the accuracy in the determination of the hamiltonian levels decreases as R and/or the energy 
increase, due to the truncation effects. The remedy to this problem consists, of course, in using a 
larger Hamiltonian, in order to reach a sufficient level of precision in the region of large R where 
we want to measure the energy splittings between the higher mass line and the threshold lines. 
In this respect, it would be interesting to make such a check by using the method employed in 
[21, where the diagonalization is performed on a larger number of states of the free fermionic 
basis corresponding to /i = 0. 
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Figure 4: First eight energy levels of the finite volume Hamiltonian of the scaling Ising model 
with magnetic field at critical temperature as functions of the scaling variable r = miR. At 
r = 40, starting from the bottom, the levels are identified as the ground state, the first three 
particle states Ai, A2 and A3, three scattering states AiAi, the particle above threshold A4. 
Crossings between the line associated to the latter and the scattering states are visible around 
r = 18, r = 25 and r = 36. 
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Figure 5: First eight energy levels of the finite volume Hamiltonian of the scaling Ising model 
with magnetic field slightly away from the critical temperature. Observe the splitting of the 
crossings pointed out in the previous figure. 
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